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We present non-linear dynamical leatures ol two-photon double-cavity optical bistability exhibited 
by a three level ladder system in the mean field limit. The system exhibits a hump like feature 
in the lower branch of the bistable response, wherein a new region of instability develops. The 
system displays a range of dynamical features varying from normal stable switching, periodic self- 
pulsing to a period-doubling route to chaos. The inclusion of two competing cooperative atom-field 
couplings leads to such rich nonlinear dynamical behavior. We provide a domain map that clearly 
delineates the various regions of stability that will aid the realization of any desired dynamics. We 
also present bifurcation diagram and the associated supporting evidence that clearly identifies the 
period-doubling route to chaos, which occurs at low input light levels. 



I. INTRODUCTION 

Understanding instability is pivotal to fabricating 
practical devices [l|. In recent times chaos has been 
used in a variety of applications including generation of 
random numbers @, Q , ultra- wide bandwidth Q , opti- 
cal communication scheme [f|, etc. Optical Bistability 
(OB) [f| 0] bas historically offered a platform for study- 
ing a variety of nonlinear dynamical effects pl-llol|. In 
conventional OB the nonlinear dynamics arises due to 
the interplay of the atom-cavity coupling as well as the 
cavity characteristics such as the number of modes, cav- 
ity decay and detuning [IJ We describe a system 
wherein the nonlinear dynamics arises due to the inter- 
play of cooperative atom-cavity coupling at two different 
frequencies. This results in a hitherto unseen nonlinear 
dynamical regime that arises in double-cavity two pho- 
ton OB, where both the fields experience independent 
feedback. This new regime occurs at low input light lev- 
els, moreover the system also exhibits negative hysteresis 
bistable response [13]. We undertake a systematic study 
that enables a greater understanding of the phase space 
structure and thus allows us to selectively steer the sys- 
tem to exhibit stable self-pulsing, chaos or regular stable 
switching. Such control mechanisms offer the possibil- 
ity of using these dynamics for communication technolo- 
gies involving multiwavelength operation, metrology, and 
moreover this system offers a different paradigm to study 
the fundamental aspects of optical instabilities. 

Optical instability has been extensively studied in the 
last four decades. Ikeda et.al. H have shown theo- 
retically that periodic instabilities and chaotic behavior 
can occur in optical bistable systems with delayed feed- 
back. Single mode instabilities ranging from gain based 
laser systems [14] to passive two-level optical bistable 
systems have been investigated earlier [lfl [l6| . Vari- 
ous studies have been undertaken related to the three- 
level atoms interacting with multiple fields leading to 
chaos: Savage et.al [l7| describe the possibility of tri- 
and quadra-stability as well as self-pulsing and chaos, 
Grangier et.al [l8| have studied OB in the purely disper- 
sive limit and shown the occurrence of chaos for large 
field intensities. Chaos has also been demonstrated in a 



three-level A system with only the probe field experienc- 
ing feedback and the coupling field detuning is used to 
drive the system to chaos [l9j . To the best our knowledge 
the regime we describe has not been reported earlier, we 
obtain chaos for moderate cooperative parameters in the 
lower branch of the bistable response at low input inten- 
sity levels. These effects crucially rely on the interplay of 
two cooperative branches within the atom. No approxi- 
mations are made with regard to the nonlinearity of the 
active medium, the effects are quite robust and occur in 
a wide parameter regime. This model is ideal for investi- 
gation of optical chaos that arises just beyond the single 
mode limit, as we consider only two single modes asso- 
ciated with the two monochromatic fields coupling the 
two-photon transition. Moreover, as the system exhibits 
coordinated dynamics (such as self-pulsing) at two dis- 
tinct optical frequencies, it can be gainfully utilized in 
optical communication technology. 

The organization of the paper is as follows: In section- 
|TT]we present the description of the theoretical model. We 
provide a detailed stability domain map as well as the de- 
tails of the numerical modelling in section- IlIII along with 
the various tools used in the analysis. In section- IIVI we 
present the results and analysis of the nonlinear dynam- 
ical features including bifurcation diagrams, spectra etc. 
We conclude in section- |Vl 



II. THEORETICAL MODEL 

We consider two optical fields in two independent uni- 
directional ring cavities that share a region containing a 
collection of three-level ladder atomic media. The op- 
tical fields couple two adjacent transitions in the atom 
with the ground state |3), and intermediate state |2) and 
the excited state |1). The electric field at the atom 
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consists of two monochromatic fields of amplitude E\ and 
i?2, at frequencies wi and L02, coupling the transitions 
|1) O 1 2) and |2) -o- |3) respectively, identical to the sys- 
tem described in the companion paper 13(. The fields 
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within the active medium are covmterpropagating in or- 
der to minimize the effects of Doppler broadening in the 
ladder system. Our analysis is within the uniform field 
approximation also known as the mean field limit. 
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FIG. 1: (Color online)The stability domain map plotted be- 
tween A2 and 62, indicating the stable fixed point region(blue) 
and the self-pulsing region(green) and the chaotic region(red), 
for the parameters 74 = 1, = 1, Ai = 0, 6\ = 0, d = 
200, I = 23, lift I = 40, for i = 1,2. The dashed line 
(62 = —8) indicates the particular parameter variation we 
consider in Fig. [2] to study the bifurcation. 

We consider the usual boundary conditions imposed 
independently on the two fields, due to the cavity feed- 
back 



atoms, di m and pi m are the dipole matrix elements and 
the atomic density matrix elements comparable to the 
|^) f-> \m) transition, and are described in detail below. 
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FIG. 2: (Color online) (a) The bifurcation diagram indicating 
the out put field |aft| and (b) The largest Lyapunov exponent 
for th same is plotted with respect to A2 for 62 = —8. The 
parameter values are same as those specified in Fig. [1] 

We undertake the analysis for the two fields along the 
lines of Ref. Q, and Eq. ^ can be reformulated using 
the above boundary conditions Eq. ([2]) as 

+ ^Pfr) (4) 

where t' = t + Atj^, the cavity decay K{ — cTi/Ci, 
the cavity detuning is scaled with «j and is given as 
Oi = 5i/ni- The absorption coefficients on along the two 
transitions are 
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for the field u>i. The subscript i = 1,2, in this paper, 
refers to the two fields at u>\ and 012, respectively. Here, 
Ti and Ri are the transmission and reflection coefficients 
associated with the two independent cavities, and L is 
the length of the active medium. The Eqs. [2] relate the 
input-output fields of the cavity. For simplicity we have 
assumed that the time taken by both fields outside the 
active medium is identical and is represented as At. The 
cavity detunings are Si — (10^ — u>i)Ci/c, where u>^ is the 
nearest resonant cavity frequency close to the field fre- 
quency u>i, and Hi is the total length of the corresponding 
cavities. The wave equation under the slowly varying en- 
velope approximation is 
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where P(cj^) is the macroscopic atomic polarization, 
which is obtained from the density matrix of the atomic 
system described below, given as P(u\) — ndi 2 pi2 and 
P(u>2) = nd 2 3P23 where n is the number density of 
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We undertake the mean field limit for both the fields 
i.e. otiL —¥ 0, Ti —¥ and Si —¥ 0, which essentially re- 
fer to spatially uniform fields that change little in each 
pass through the cavity, however, owing to the large cav- 
ity photon lifetime the photons undertake many passes 
inside the cavity [ill ] resulting in significant cooperative 
effects. The governing equations of motion involving the 
scaled input and output fields are 



dx 1 
dx 2 



T = K 1 [-x 1 {l + i6 1 )+y 1 + 2iCipi2] , (5) 
T = K2[-x 2 (l + ie 2 ) + y2 + 2iC2P23} ■ (6) 



where the dimensionless input /output fields are de- 
fined as 
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FIG. 3: (Color online) The bifurcation points associated with 
the bistable response for the field at u;2 while the input field 
at uji is held constant. The parameters are C\ = 200, C2 = 
200, = 23. Inset: The bistable response in absence of 
feedback to field at u\ (Ci = 0, C2 = 200) with \y\ \ = 1.5 and 
other parameters are 7* = 1,/tj = l,Ai = 0, 9\ = 0, A2 = 
4,02 = -3 where i = 1,2. 



The last terms in Eqs. (|5I6[) contain the atom-cavity 
coupling via the cooperative parameters C\ — a\L/2T\ 
and C2 = OL1LIIT2 along the two transitions. The atom- 
field interaction is governed by the density matrix equa- 
tions of motion 
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-iG* 2 Piz 
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— 1G2P12 
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-(72 + ?A 2 )p23 + iG 2 {p33 - P22) 

+iG* lPl3 

272P22 - iG 2 p32 + iG* 2 P23 (10) 



where the atomic detunings are given by Ai = wi 2 — ^1 
and A 2 = UJ23 ~^2-, and 271 and 272 are the spontaneous 
emission rates from the levels |1) to |2) and |2) to |3), re- 
spectively, thus considering only the radiative relaxation 
processes in the dilute atomic gas active medium. The 
schematic of the level scheme as well as the the two cav- 
ities are given in the companion paper All the fre- 
quency units are normalized with respect to the atomic 
decay 72, unless specified otherwise. Furthermore, we 
undertake a detailed non-linear dynamical analysis of the 
system described above. 



where the total Hamiltonian H is 



III. NUMERICAL MODELLING AND 
STABILITY MAP 



H — H a t + Hint , 

H at = fiwi 3 |l)<l|+fej23|2)<2|, (9) 
H int = -d-E = -hGx\l)(2\-hG2\2)(3\+h.c. 

The terms Hu>i3(hui23) correspond to the energy levels 
of bare atom measured from the ground state |3), and 
the interaction Hamiltonian Hint is given in the dipole 
approximation, involving the Rabi frequencies 
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The relaxation processes like spontaneous emission and 
dephasing of the atomic coherence are contained in the 
Liouville operator C of Eq. (jHJ.We explicitly enumerate 
the equations of motion of density matrix elements under 
the rotating wave approximation 



In order to understand the nonlinear dynamical as- 
pects of this system we solve the field equations Eqs (|5|6j) 
self consistently with the atomic evolution given in 
Eqs. (jlOp . Investigation of the non-linear dynamics is car- 
ried out by computing the response (such as ouput fields 
Xi) for a given set of input fields y\ and j/2, for a certain 
choice of the atomic parameters such as A; and 7j, the 
cavity parameters 9i, and Ki and the atom-cavity cou- 
pling governed by the cooperative parameters CV The 
Newton-Raphson method is used to obtain the solutions 
of the coupled nonlinear system of equations in the steady 
state limit. The atom-field cavity interaction dictates the 
phases and amplitudes of both the output fields. The 
multiplicity of the underlying solutions needs a careful 
handling with regard to its numerical computation and 
is described in detail in the companion paper [l3j . In 
essence, this system does not permit an apriori choice 
of the output fields, which is undertaken conventionally 
to compute the corresponding input field (in systems in- 
volving feedback for one field). 

The computations were undertaken using two indepen- 
dent numerical tools, one involving Fortran libraries asso- 
ciated with EISPACK [H and the other using MATLAB. 



4 




liol 


Period-one 


Period-two 


Period-four 


49.78 
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-1.77 + 0i 
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TABLE I: The Floquet multipliers for the first few points of 
periodic doubling cascade for the input field indicated in the 
inset. (ii) of Fig. Q|a) 
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FIG. 4: (Color online) Bifurcation diagrams of the cavity out- 
put fields, (a) \x2\ versus |j/a| and (b) versus \y2\, for the 
same parameters as in Fig. [3] 



Furthermore, we have also used the MATLAB continua- 
tion package MATCONT j2l| for study of the nonlinear 
dynamical aspects of this system. We identify a series of 
fixed points that exhibit interesting bifurcations such as 
Hopf points, Limit points as shown in Fig. [3l 

We present the highlights of the linear stability analy- 
sis in the form of a detailed stability map shown in Fig.ffl 
The unstable regions and its neighborhood are quite often 
related to a variety of dynamical behaviour such as peri- 
odic, quasi-periodic or even chaotic dynamics, and hence 
a detailed bifurcation analysis is undertaken. Specifically, 
the periodic self-pulsing dynamics is further character- 
ized by the Floquet multipliers that distinguish between 
the stable and unstable limit cycles, and the Lyapunov 
exponents are calculated in order to identify and charac- 
terize chaotic behavior [l[ . All these nonlinear dynamical 
features are independently confirmed by integrating the 
time-dependent equations Eqs. ©,([6]), and ( [10]) using 
the fourth order Runge-Kutta technique p2j |. 

We believe the stability domain map would facilitate 
the experimental realization of these diverse nonlinear 
dynamical features. The parameter space associated with 
the system is exceedingly large owing to twelve differ- 
ent physical parameters 71,2, ^1,2, 61,2, «a,2i Ci,2, E\ n 2 
each of which can be varied independently. The stability 
map in Fig. Q] indicates the range of just two of these pa- 
rameters A2 and 6*2 for which the system exhibits stable 
switching (blue region) , periodic or quasi-periodic (green 
region) and chaotic (red region) . One can clearly see the 
islands (green) of stable self-pulsing within the chaotic 
region (red). The dashed black line indicates the dy- 
namics that we further explore and is discussed below 
in detail along with the associated bifurcation diagram. 
The stability domain map allows one to pickup the ap- 



propriate parameter regime that corresponds to a desired 
dynamical behavior, such as self-pulsing, stable switching 
or chaotic output, even a complete avoidance of chaotic 
region is possible by a judicious choice of the trajectory 
in the multidimensional parameter space. 

IV. NON-LINEAR DYNAMICS 

In order to understand the route to chaos we consider 
the bifurcation diagram along the dashed line in Fig. [T] 
which is shown in Fig. [2ja). As one varies the detuning 
of the field A2 while maintaining a constant cavity de- 
tuning, we obtain a period doubling cascade for the val- 
ues of A2 within [0 : 1.2] and inverse periodic doubling 
around A 2 ~ 8.5. The largest among the Lyapunov ex- 
ponents corresponding to the regions of self-pulsing are 
close to zero. The chaotic domains are clearly identified 
by the positive largest Lyapunov exponents as shown in 
Fig. EJJb). The self-pulsing and chaotic dynamics occur 
intermittently and the windows of stable periodic self- 
pulsing within a largely chaotic region are clearly seen 
in Fig. [2jb) . Beyond A2 ~ 8.5 all the Lyapunov coeffi- 
cients become negative indicating that the system quickly 
moves towards a stable fixed point behavior independent 
of the initial conditions, in our case it corresponds to 
stable switching. 

We illustrate the essential nonlinear dynamical fea- 
tures where the magnitude of the input field yi = 23, 
while the other incident field 2/2 is varied. In the S- 
shaped bistable response we obtain a hump like feature, 
and the associated linear stability analysis is shown in 
Fig. [21 black (red) color indicate unstable (stable) fixed 
points. The corresponding bistable curve without the 
cavity feedback for the field coupling upper transition 
(systems with cooperative parameter C\ = 0) is shown 
in the inset of Fig. [3] This hump like feature under 
appropriate parameter regime transforms into negative 
hysteretic bistable response. Furthermore, this regime 
is also associated with non-linear dynamical behavior, 
and is intimately related to the cavity assisted inversion 
as discussed in the companion paper [l~3l ]. The hump 
exhibits Hopf bifurcations and these effects are absent 
in OB models involving feedback for one field(inset of 
Fig. however we note that, this feature is a result 
of an intricate interplay of phases of both the fields. If 
one demands that both the output fields are real, as in 
Ref. [23j], this hump like feature as well as the nonlinear 
dynamical features disappear. We indicate the bifurca- 
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FIG. 5: (Color onlone)Time evolution of output field X2 at uu2 are shown in the first column, the corresponding spectrum 
(in arbitrary units) and the population phase plots (pn versus P22 and P33) are shown in the middle and the last column 
respectively. Plots in each row corresponds to various values of 1 3/2 1 involved in the sequence of period doubling route ultimately 
leading to chaos, (a) Period one with I2/2I = 49.78, (b) Period two with \y2\ = 49.01, (c) Period four with 1 2/2 1 = 47.70, and (d) 
Chaos with | j/2 1 = 32.47. The other parameters are same as those for Fig. [3] 



tion points such as the Hopf points(H), limit points(LP) 
in the figure. In the unstable regime corresponds lower 
cooperative branch, one obtains self pulsing as well as 
chaos at low input light levels. The bifurcation diagram 
associated with Fig. [3] is shown in Fig. @] We also ex- 
pand two illustrative regions of the bifurcation diagram 
that indicate the onset of chaos. With increasing input 
field 2/2 the system loses stability of stable fixed point so- 
lutions at first Hopf point (Hi) and exhibits self-pulsing 
behavior. This is indicative of a supercritical Hopf bi- 
furcation at Hi. The Hopf points indicate the onset of 
periodic behavior, and in between the Hopf points we ob- 
serve a period doubling route to chaos. The stability of 
the periodic behavior is established using Floquct multi- 
pliers, and chaos is confirmed using Lyapunov Exponents 
as well as the output field spectrum. 

There is a simultaneous existence of stable fixed point 
solutions as well as self-pulsing limit cycles beyond the 
second Hopf point Hi- For certain range of yi magni- 
tude within (48 : 55) depending on the initial condition 
the system can be driven to exhibit either periodic self- 
pulsing or regular stable switching. A transition from 
periodic behavior to fixed point solution behavior occurs 
as a Saddle- node/fold bifurcation takes place beyond H2 
a t J/2 = 54.6, and is represented as limit point of cycles 
(LPC) (see Fig. g]). At LPC two limit cycles (one sta- 
ble and one unstable) coalesce and annihilate each other, 
leading to pure fixed point solutions [l[ . The stability of 
limit cycles indicated by the Floquet multipliers remain- 
ing with in unit circle and transition out of unit circle 
implies unstable limit cycles. The Floquet multipliers 



associated with the stable and unstable limit cycles (at 
1 2/2 ] = 51 before coalesce at LPC) are 0.3639 + Oi and 
— 1.7997+ Oi respectively. 

The onset of chaos indicated in Fig. U] occurs as one 
decreases the input field j/2 magnitude from 50 (see in- 
set (ii)). At first one obtains self pulsing output of 
period 1, which transforms into a self pulsing output 
of period 2 and then to period 4, as shown in Fig. [S] 
(a-c). The genesis of new frequencies in the spectrum 
of the output field (second column) at fp2 = /pi/2, 
/p4 = /pi/4 and its multiples, apart from the dominant 
frequency fpi — 0.48, clearly indicates a period doubling 
route. The dashed arrows in Fig. indicate conventional 
bistable switching. As both the fields switch simultane- 
ously, one can utilize this system for controllable switch- 
ing at two distinct optical frequencies in optical com- 
munication applications. The period doubling route to 
chaos is further substantiated by the three dimensional 
phase plots (third column) involving the population of 
the three atomic states which show closed curves involv- 
ing one, two and four loops, respectively. As we further 
decrease the input field 7/2 magnitude to 32.47 (indicated 
by thick line with cyan color in Fig. EJa)) one observes 
chaos and the corresponding two largest Lyapunov expo- 
nents of the system are 0.3356 and 0.0031. The spectrum 
of the output field becomes continuous as well as the pop- 
ulation phase plot gets filled up clearly establishing the 
existence of chaos. This is further corroborated by the 
loss of stability of periodic behavior (limit cycles) [24| as 
shown in Table. HI wherein the associated Floquet mul- 
tipliers cross out of the unit circle along the negative 
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real axis. The Floquet multipliers having zero imaginary 
part and transition out of the unit circle along negative 
imaginary axis precludes the existence of quasi periodic- 
ity, thus establishing a periodic-doubling route to chaos 
as indicated. The nonlinear dynamics of both the fields 
closely resembles each other, and one obtains periodic 
self-pulsing to chaotic dynamics for both the fields X\ 
and X2 simultaneously. 




FIG. 6: (Color online) The development of negative hysteresis 
and multistability for the field at 0J2 with Ai = 4, A2 = 0. 
Each curve corresponds to a particular value of the input field 
at uj\ whose magnitude is varied from 20 to 30, and the other 
parameters are same as in Fig. [3] 

This system offers a variety of control mechanisms 
which allows one to access any desired dynamics by 
merely changing a few pertinent parameters. The in- 
clusion of finite atomic detuning Ai could even lead to 
multistability as seen in Fig. [5] The figure also indicates 
the various domains of stability where interesting non- 
linear dynamical behavior occurs. A detailed analysis of 
the multistable behavior will be presented elsewhere. 
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FIG. 7: (Color online) (a)The bifurcation diagram for the 
field at LJ2 as the ft 2 cavity decay is tuned. (b),(c) and 
(d) shows the time evolution of output field magnitude for 
k 2 = 0.7808, 0.3008, 0.5803 (marked as i, Hi and ii in (a)) for 
C x = 1000, C 2 = 100, = 25, \y 2 \ = 13, 2 = 0, A 2 = and 
the other parameters are same as in Fig. Q] 

We indicate other ways to control the stability as well 
as nonlinear dynamical aspects in this system. One can 
obtain negative hysteretic behavior, with the associated 
non-linear dynamical features in its neighborhood. One 



can obtain all the above mentioned non-linear dynami- 
cal features by varying the cavity decays «i and/or K2- 
Such control over the nonlinear dynamical response is 
illustrated in the bifurcation diagram shown in Fig. [7] 
which involves K2 as the control parameter which can be 
used to switch on or off the chaotic dynamics. It should 
be noted that including finite cavity detuning (#2 ) aids 
in obtaining all the non linear dynamical features with 
comparable values of the co-operative parameters C\ and 
C2 as shown in Fig. [H 




FIG. 8: (Color online) The hump like feature progressively 
transforms into a negative hysteresis for comparable values of 
Ci and C2, and other parameters are same as Fig. [3] 



V. CONCLUSIONS AND REMARKS 

We have demonstrated a new regime of nonlinear dy- 
namical response at low input light levels for two-photon 
double cavity optical bistability with three-level ladder 
atomic system as the active medium. Independent feed- 
back is applied for both the fields interacting with the 
atom. We present bifurcation diagrams that allows one 
to access the desired regime of dynamics for both the cav- 
ity fields. We prove the period-doubling route to chaos 
in this new regime associated with the lower cooperative 
branch. A control paradigm based on careful maneu- 
vering of parameters so as to traverse across phase space 
in order to obtain any desired dynamics is demonstrated. 
The system exhibits negative as well as positive hysteresis 
bistable response, stable periodic self-pulsing and chaotic 
dynamics, apart from the conventional normal switching 
and multistability. 
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